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In this paper we show that both music composition and brain function, as revealed by the Elec- 
troencephalogram (EEG) analysis, are renewal non-Poisson processes living in the non-ergodic do- 
minion. To reach this important conclusion we process the data with the minimum spanning tree 
method, so as to detect significant events, thereby building a sequence of times, which is the time 
series to analyze. Then we show that in both cases, EEG and music composition, these signifi- 
cant events are the signature of a non-Poisson renewal process. This conclusion is reached using 
a technique of statistical analysis recently developed by our group, the Aging Experiment (AE). 
First, we find that in both cases the distances between two consecutive events are described by 
non-exponential histograms, thereby proving the non-Poisson nature of these processes. The corre- 
sponding survival probabilities ^(t) are well fitted by stretched exponentials cx exp(— (yt) a ), 
with 0.5 < a < 1.) The second step rests on the adoption of AE, which shows that these are renewal 
processes. We show that the stretched exponential, due to its renewal character, is the emerging 
tip of an iceberg, whose underwater part has slow tails with an inverse power law structure with 
power index fi = 1 + a. Adopting the AE procedure we find that both EEG and music composition 
yield /i < 2. On the basis of the recently discovered complexity matching effect, according to which 
a complex system S with /is < 2 responds only to a complex driving signal P with fip < fis, we 
conclude that the results of our analysis may explain the influence of music on the human brain. 

PACS numbers: 05.40.-a, 02.50.-r, 87. 19. La, 89.75.Hc 



I. INTRODUCTION 



The study of neuronal systems is a challenge for statistical physics, insofar as experimental evidence is proving 
that the ordinary Poisson paradigm is inadequate to deal with these complex systems. According to some neuro- 
physiologists the neurons are renewal [l[ and they are markedly non-Poisson More precisely, the experimental 
evidence of in vitro observations coupled with analysis of in vivo spiking patterns indicate that single neurons are 
fundamentally non-Poisson processes In the literature the implicit assumption is frequently made that even if 
the spiking activity of a single neuron is not Poissonian, the activity of a set of many neurons is Poissonian. This 
assumption may lead us to conclude that the human brain is a Poisson system. The authors of the paper of Ref. [H 
proved that this assumption is invalid. In this paper, through the Electroencephalogram (EEG) analysis we reach 
the conclusion that the human brain is not a Poisson system, in line with the theoretical remarks of Ref. pj. It 
is important to point out that in this paper the term non-Poisson process indicates a strong deviation from the 
exponential decay. This kind of non-Poisson behavior implies, as we shall see, the emergence of fat tails with an 
inverse power law behavior, although for many reasons, ranging from the finite time observation to the influence of 
spurious random fluctuations, in the long-time region these tails are truncated. 

In this paper we apply the same statistical analysis to music composition. There is a wide agreement that music 
composition can be thought of as a complex signal. It is convenient to quote the seminal work by Voss and Clarke [f|. 
These authors have found in fact that music composition yields a 1// noise spectrum, which is generally regarded as 
a complexity manifestation [f| 0, H, [§] . 

The power of music to evoke emotions is well known, but only recently, as discussed in the short review of Ref. 
[Tfj| . it has attracted the attention of neuro-scientists [i ll. How does the communication between music and the 
brain take place? The authors of an interesting paper [12| have recently studied the brain response of musicians and 
non-musicians to music listening and have found that musicians yield a higher degree of the gamma band synchrony. 
A more recent study of these investigators [l3j has supported the theory that phase synchronization is a significant 
marker in human cognition (l4| . 

The phase synchronization established by the authors of Refs. [HI, [l3[ is one of the interesting properties of chaos 
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synchronization [T^|, which is in fact attracting the increasing interest of neuro-scientists According to the 

seminal work of Ref. [l7j both driving system (the music com pos ition) and response system (the brain) are dynamic 
systems with strange attractors (for instance Rossler systems [fa]). 

In this paper we analyze both EEG data and music composition with the joint use of two techniques, the Minimum 
Spanning Tree (MST) approach and the Aging Experiment (AE). Both EEG and music composition data are expressed 
as multidimensional vectors Vi, with i = l,...,N. The different components of these vectors correspond to the 
information afforded by different electrodes in the EEG case, and timber, pitch, harmony, melody, rhythm . . .etc., 
in the case of music composition. In this representation the existence of events is not evident. Thus to reveal their 
presence we build the MST proposed by Kruskal and we study its time evolution: From time to time the topology 
of this MST undergoes an abrupt change that we interpret as an event. The probability of an event occurrence at a 
time distance r from the preceding event is found to be a stretched exponential. With the help of the AE method 
we prove this process to be renewal. The conclusion of this paper is that both EEG and music composition are 
non-Poisson renewal processes, generated by an inverse power law distribution with index fj, < 2. This condition 
apparently departs from those conditions considered to be essential for synchronization, and so for the brain to be 
sensitive to music. We note however that it has been shown [13, H3 that while a non-Poisson complex system does 



not respond to harmonic perturbation \2'2\ . 23 1 , it is very sensitive to the influence of a complex perturbation with 
the same complexity. The authors of Ref. 2C, [2I| denoted this effect as Complexity Matching (CM) phenomenon. 
Therefore, the non-Poisson renewal condition is compatible with the transmission of information from a complex 
system to another, and, on top of that, thanks to the CM effect, might allow us to explain the influence of music on 
the brain. 

The outline of this paper is as follows. In Section [TlJ with the help of Appendix A, we show under which condition 
a survival probability (SP) with the form of a stretched exponential can be considered a non-Poisson renewal process 
with power law index fj, < 2. In Section HTTl we illustrate the MST method. Sections llVl and fVl illustrate how we reached 
the conclusion that both EEG and music composition data are manifestations of non-Poisson renewal processes with 
fi < 2. In Section IVll we concisely review the CM effect so as to support the main conclusion, illustrated in Section 
IVII1 that it is possible to tune music composition to the brain so as to make the brain respond to music. 



II. THE STRETCHED EXPONENTIAL AS A TRUNCATED MITTAG-LEFFLER FUNCTION 

The purpose of this Section is to derive the following conclusions: (i) A particular class of stochastic processes 
leads, in an appropriate limit, to a renewal process with SP described by a Mittag-LefHcr (ML) function; (ii) The ML 
function has the small r limiting form of a stretched exponential, exp(— (71"")), with a < 1, and a large r power law 
limit; (iii) The power law index is related to the stretched exponential parameter a by Eq. (iv) One can deduce 
the power law index from the small r stretched exponential fit, even in situations where there is no information about 
the long tail in the data recorded. This is important for the main purpose of this paper. In fact, in Sections IIVI and fVl 
we shall prove that EEG's and complex sounds, respectively, are examples of such a process because they are renewal 
processes and are well fitted by a stretched exponential. By means of Eq. ^ we shall derive the hidden information 
on fi. 

The Mittag-Leffler (ML) function [24| is attracting an ever increasing interest in the literature of complex fluids 
such as liquid crystals, glass- forming liquids, and polymeric and colloid systems [25| . A remarkable property of the ML 
function is that the complex susceptibility produced by the ML relaxation function yields the Cole-Cole experimental 
form 26]. For this interesting property the reader can consult also Refs. [25|, [13] and [28| . 

The ML function establishes a bridge between the stretched exponential behavior for short time and an inverse 
power law in the long-time limit [29j | . This property is important, since it may settle the controversy between the 
advocates of stretched exponentials and those of inverse power laws. In the case of the financial market, the authors 
of Ref. [3(| have found that the ML function affords a very good fitting of experimental data. There is a problem 
with the tails, however, insofar as the lack of sufficiently rich statistics make noisy the time region where the ML fat 
tail should appear [3I( . The authors of Ref. [3I( used the aging experiment to make an inverse power law behavior 
distinctly emerge from the noisy background of the long time regime. 

Here we show an approach to the ML proving that in some noisy conditions only the stretched exponential portion 
of the ML function can remain visible. Our theoretical approach refers to the time distance between two consecutive 
events, rather than a molecular relaxation process (26|: The motor driving our process is a physical generator of 
events. Let us assume that a natural time scale exists, where this physical generator of events is a Poisson process 
with rate r <C 1, so as to make the discrete time representation virtually indistinguishable from the continuous time 
picture. The time distance between two consecutive events produced by this generator is given by the waiting time 
distribution density ipp(n), which has the following exponential form 



4>p(n) — rexp(~rn). 



(1) 
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To generate the sequence of events of interest for this paper, we operate as follows. We record the activity of this 
physical generator at each and every time step of its natural time n. At the same time, following the prescription of 
the subordination theory (3^ |. we turn the time n into the continuous time t(n), by setting 

t(n+l)-t(n)=T n . (2) 

It is straightforward to prove (see Appendix [Aj that the adoption of an inverse power law with no truncation and 
power index [i as subordination function, together with the renewal assumption for the events production and the 
condition r«l, creates a SP with the form of a ML function of order a, denoted by the symbol f jp, whose analytical 
expression is 

with 

a = n-l. (4) 
In the case < a < 1 it is well known that the ML function admits two limiting conditions, namely 

£ a (-( 7 *n~^ t-oo, t> 1 - (5) 

that is an inverse power law, and 

E a (-(lt) a ) ~ exp(-(jt) a ) i-0, t<-, (6) 

7 

which is a stretched exponential. 

To take into account that the data under study are finite, it is convenient [32j to select for the subordination 
function ip(r) an inverse power law with index /i, which is exponentially truncated at t > where L < 7. As an 
effect of this choice (see Appendix [AT) we get for the function ^sp(t) a resulting form virtually indistinguishable from 
a stretched exponential, in the intermediate time scale where the departure from the exponential form is significant. 
This derivation of the stretched exponential is different from the one recently proposed by the authors of Ref. [33j ] . 
However, a comparison between our approach and that of these authors is not quite appropriate: in fact, Magdziarz 
and Weron [33[ aim at the same purpose as that of the earlier work of Ref. (26|, the explanation of the non- 
exponential Cole-Cole relaxation, and do not afford prescriptions to evaluate the distribution of the time distances 
among consecutive renewal events, which is the main trust of this Section. 

The subordination process is realized by means of a random prescription and, consequently, ^sp{t) is expected to 
fit the renewal condition. We shall assess this property by means of the aging experiment of Section IIVI We shall 
also establish a connection between the stretched exponential exp(— (jt) a ), with a < 1 and the power index \x by 
means of the relation of Eq. ^ generated by the subordination procedure. Since the subordination is realized with 
a function ^(t), which is a truncated inverse power law, the resulting process in the long-time scale does not violate 
the ergodic condition. Thus the Poisson behavior, and the consequent lack of aging, as we shall see by means of the 
statistical analysis of both EEG data and music composition, is recovered. The same limitations are shared by many 
other complex processes, see for example [32], due to the obvious fact that an exact inverse power law behavior is an 
idealization that would imply the infinite size of the systems under study. The study of this idealized condition is 
useful to shed light into the transient behavior before the eventual Poisson condition. 



III. THE MINIMUM SPANNING TREE APPROACH AS A GENERATOR OF EVENTS 



In this Section we introduce an algorithmic procedure that will allow us to process the data at our disposal so as to 
generate events and consequently time series to analyze. The method is based on the famous Minimal Spanning Tree 
(MST) algorithm [34j | , which we now briefly introduce. 

To define the MST approach we closely follow the arguments of Ref. [Hj]. Imagine a data set consisting of N 
columns, each column representing the signal recorded by an electrode, or the timber, pitch, note . . . , etc, of the 
music composition. Taking two columns, say x and y, we define the correlation coefficient between the two columns 
as follows 

„ _ SLofofc- < x >)(vk- <y>) f7 , 

Pxy — , { I ) 

X Oy 
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where the quantities < x > and a x are, respectively, the average and the standard deviation of the values that x 
takes over the interval [0, t]. Consider now the following quantity 

x— < x > 

x = (8) 

Ox 

and assume that in the interval [0, t] there are t values of the vector x, namely £k, k — 1, . . . , t. The distance between 
two columns x and y over the time interval t is easily obtained adopting the well known formula of the Euclidean 
distance 

t 

dl y = \\x-y\\ = J2^-y^- W 



k=l 



Moreover J2k=i ^fe = 1- Therefore, Eq. (J9j) becomes 



t 



dl y =J2(il+yt-^ k yk) = 2-2j2ikyk- (10) 



k=l k=l 



The last sum on the right hand side of Eq. (JTDJ) is the correlation coefficient between x and y of Eq. ((7} at the time 
t. It is therefore possible to define the distance between the columns x and y as 



d X y=\j2{l- Pxv ). (11) 

Since — 1 < p xy < 1, then < d xy < 2. Moreover d xy fulfills the property of a distance [35|. 

Through the distance d xy we now introduce the MST approach over the time interval t. In a connected graph of 
N objects, weighted through the distances d xy , the MST is the tree with N — 1 links for which the total sum of 
the edges is minimum. In the literature there are many algorithms to create the MST. Here we select the method 
proposed by Kruskal [191 ]. This method consists of the following steps: we sort the distances in increasing order; we 
select the shortest distance and draw an edge between the associated nodes; we go to the next distance and draw 
edges; if an edge creates a loop we erase it; we continue drawing lines until all the N columns are represented. 

In this paper we make a dynamical use of the MST. We build the MST over a time windows of length t, then 
we move to the next (non overlapped) window and build the MST corresponding to the new position. An event 
is defined whenever a change in the distribution of links occurs. This way of proceeding allows us to define events 
and consequently the distribution of waiting times between two consecutive events. Moreover, when building the 
distribution of the number of edges, we do not label our electrodes, but simply look at the shape of the distribution, 
that essentially characterizes the topology of the MST. If distributions at subsequent times do not exactly coin- 
cide, we consider it to be a signal of a change in the global properties of the brain, and indicate this as our critical event. 



The MST has been recentl y a dopted in the literature in an increasing number of papers on different complex 
systems, from networks (e.g. [36j | among the others) and financial markets (e.g. [37j among the others) to neuro- 
physiological processes [38[ . Particularly, we quote the recent paper of McDonald et al. [39| . In this paper the MST 
approach has already been used and applied to financial time series, but in a different way. In fact, the authors 
of Ref. monitor the survival-ratio of edges in time, that is, the ratio between the number of edges connected 
to a certain node at time and the same quantity computed at a later time t, and find that it decays non exponentially. 

According to Ref. [3^|, the MST, by definition, is more suitable to study positively than negatively correlated sys- 
tems. Indeed, negatively correlated electrodes contribute differently to the MST than positively correlated electrodes. 
In our analysis this criticism only relatively applies, as we arc interested in the changes of the MST topology in time, 
i. e. in the changes of the correlations with time. 

In the next Section we shall analyze the statistical properties of the time series obtained in the aforementioned way. 



IV. ANALYSIS OF EEG 



Digital EEG data was collected on 5 subjects, utilizing a Scan LT-40 amplifier, manufactured by Neuro Scan 
Medical Systems, El Paso, Texas. The Scan LT-40 is an FDA approved medical device for the collection of digital 
EEG. Online monitoring of EEG was provided by NeuroScan Medical Systems 1.2 software. Offline evaluation and 
the removal of artifacts from the EEG record and the conversion of the data to ASCII files were accomplished using 
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NeuroGuide 2.2.6 software. An electrode cap from Electrocap International Inc. was used to provide standardized 
electrode placement. Digital EEG data was collected from 19 locations using the international 10-20 system 
of electrode placement. A reference electrode was placed on each earlobe to provide a linked ears montage for 
the physical reference of the scalp recordings. The impedance of the respective earlobe reference electrodes was 
maintained within 1 KVl of each other. All other electrode impedances were maintained at 1 or 2 KQ, relative 
to amplifier input impedance with no more than 1 Ktt of variance between any of the electrode contacts. The 
amplifiers used to acquire the EEG were calibrated with sine waves before the acquisition of EEG. The EEG data 
was digitized at a rate of 250 samples per second. Before Analog to Digital conversion, anti-aliasing was achieved 
by a low-pass filter built into the software. The EEG was visually inspected online during acquisition to monitor 
for artifact. When necessary, data collection was stopped to identify and remove persistent sources of artifact such 
as muscle tension. The subjects included 2 healthy subjects, 1 subject with chronic back pain, and 2 subjects with 
mental depression. Subjects were medication free during the data acquisition. The average number of data collected 
is 26, 141, corresponding to about 104 seconds of record. The minimum number of data for a single subject is 10, 961, 
corresponding to about 43 seconds of recording, for one of the healthy subjects, while the maximum number of data 
for a single subject is 56, 544, corresponding to about 226 seconds, for the subject with chronic back pain. Only two 
states were considered for the subjects, namely eyes open (EO) and eyes closed (EC). The acquisitions were made in 
the same conditions for all the subjects. 

The MST approach, discussed in the previous Section, is applied to the set of 19 columns that the data acquisition 
method affords. The time series of MST topological changes is obtained and, to be consistent with the theoretical 
results of Section|ITl the SP of the time distances between two consecutive MST topological changes is evaluated. 
In Fig. [T] we show the SP relative to the individual with back pain in EO condition. We fit the resulting curve with a 
stretched exponential of the form of Eq. © with parameters 7 = 0.205 and a = 0.595. 

1 
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FIG. 1: The SP of MST topological changes for the individual with back pain in the EO condition. The curve is fitted using a 
stretched exponential of the form of Eq. ([6]) with 7 = 0.205 and a = 0.595. 

Similar results have been obtained for the SPs of all the other individuals at our disposal. Using Eq. (j4|), we find 
that the values of /j, for the EO condition lie between 1.595 ± 0.005 and 1.960 ± 0.013, while for the EC condition 
between 1.746 ± 0.018 and 1.974 ± 0.024. 

The results of Fig. [T] suggest that the process is strongly non-Poissonian. However, the mere analysis of the 
distribution of waiting times is not enough to establish the real nature of the process. In fact, as pointed out in 
Refs. [13,|4l[, a modulated Poisson process as well can have as an outcome a non-exponential waiting time distribution. 
In order to assess if the time series under study is produced by a genuinely non-Poisson renewal process, the authors 
of Ref. [4l[ propose the Renewal Aging Experiment (AE). According to this procedure, in addition to the waiting 
time density ip(r), an aged waiting time density il)exp(i~,t a ) must be evaluated. This is done as follows. Let {U} be 
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the series of time obtained through the prescription indicated in the previous Section. For each time tj the first time 
of the sequence at a distance from ij equal to or larger than ij + t a is recorded. This time will be t^, with k > i. The 
time distance r(ti,t a ) — t^ — (i, + t a ) is considered. The procedure is repeated for all the times of the sequence {U}, 
and the sequence of these recorded time distances is used to generate the distribution density ip exp {T, t a ). Moreover, 
the following quantity is evaluated 



Jo" dyip num (T + y) 
K(t a ) : 



(12) 



where K(t a ) is a suitable normalization constant and ip nwn (t) is the numerical waiting time distribution corresponding 
to t a = 0. 

To establish the renewal character of the process by means of the AE it is convenient to plot the corresponding 
SP's, ty ren (T, t a ) and ^> exp (T, t a ). The SP of age t a — 0, indicated as \&o, is also plotted. In Fig. [2] the three SPs are 
compared, and it is found that ^ e x P and $ re „ virtually coincide. According to Ref. [4l| this good accordance is the 
numerical evidence of the renewal nature of the process. 
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FIG. 2: The AE analysis on the time series of Fig. [T] for t a = 30. The aging is present and compatible with the renewal 
condition. 



Moreover this result confirms the theoretical framework of Section lU allowing us to obtain the index of the inverse 
power law distribution from the index of the stretched exponential, namely from Eq. ^ . Since all the values of a are 
smaller than 1 we conclude that /i < 2, and therefore that the human brain is a non-ergodic system. It is important to 
repeat here what we already pointed out in Section II, that in the long run, as all the anomalous processes occurring 
in nature, the ergodic properties, as well as the Poisson condition and the associated lack of aging, are recovered. 

How can we explain the emergence of this stretched exponential of renewal origin? An unknown referee pointed 
to us that there may exist a dominant structure and that the spanning trees appearing in the processed data may 
produce short and long waiting times according to their distance from the dominant structure. If there is a systematic 
alternation of short and long times, we end up in a model recently discussed in Ref. |40| . In his case there would 
be a drastic reduction of renewal aging. If on the contrary the process is not a systematic alternation of short and 
long waiting times, it becomes indistinguishable from the non-Poisson model that we propose with the theory of 
Section HU 



Now we have to address an important issue, concerning the physical origin of the complexity property revealed by 
our analysis. We share the opinion of Haken [42| that the brain global behavior is an emergent property produced by 
a synergetic process of synchronization [43| ■ As a consequence the different electrodes are correlated and the MST 
of Section III detects the emergence of this global property. The electrodes on which the MST analysis is based, 
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are single units, like the interacting columns of a surface growing as an effect of random deposition [32J. Due to the 
cooperation with the other columns the single columns inherit the complexity of the whole growing surface [32| . 

On the basis of this observation, we are tempted to make the conjecture that the single electrodes inherit the global 
complexity. It is important to stress that Buiatti and his co-workers |44| have recently obtained results confirming 
this conjecture, with a method of analysis based on the observation of a single EEC Our conjecture is also compatible 
with the results found by the authors of Ref. (4f| . To understand the connection between the results of this section 
and those found by the authors of Ref. [45| the reader should consult the recent article of Ref. [46j . The authors of 
this paper have studied the dynamic approach to Fractional Brownian Motion (FBM), as expressed by 

j t x{t)=m, as) 

where is a fluctuating velocity with memory, namely, a correlation function with slow tails. The dynamical 
approach to FBM proves that the variable x in the asymptotic time regime shares the same properties as the traditional 
form of FBM. They have also proved that the origin recrossing of x(t) generates a non-Poisson renewal process with 
the power /i related to the FBM scaling coefficient H by 

fi = 2-H. (14) 

Note that the adoption of the more realistic model of the interacting columns of an interface growing as an effect of 
random deposition of particles would produce saturation and also a truncation of the inverse power law waiting time 
distribution, a property similar to the long-time Poisson behavior revealed by Fig. 1 and Fig. 2. However, if the 
observation is limited to the short-time region, the condition of Eq. (Ti"4"|) applies. 

The analysis of the authors of Ref. [45| rests on £(t), rather than on x(t). We are convinced in fact that the single 
electrodes yield a signal that has to be interpreted as the variable x of Eq. (fT3|) . Thus to relate the results of Ref. [45| 
to the results of this Section, we must ado pt E q. (TT4")) . This procedure generates fi < 2. With this interpretation in 
mind, we find that also the analysis of Ref. [45[ proves that individuals whose signal was recorded in the EC condition 
is closer to /i = 2 than the individuals in the EO condition. 

Stressing this result is a way for us to draw the attention of the reader on the fact that the discovery of non-Poisson 
renewal events, made possible by the method adopted in this paper, does not yield results conflicting with the work 
of other groups. In Table HI we present a sample of the results of our analysis on our group of individuals. We remind 
that two of them were healthy, one had back pain, and two had mental depression. The values of [L refer to the same 
subject, in the two conditions. We see that the parameters y, of individuals in the EC condition are significantly closer 
to 2 than the parameters fi of patients in the EO condition. In Table |TT] we report a sample of the results of Ref. 
PH (healthy individuals only, Tab. I) in terms of fi, by means of the rule of Eq. (|14p . By comparing the results of 
Table U with those of Table |TT] we reach the conclusion that the results of our analysis are compatible with those of 
Ref. [43]. It is necessary to stress that this result is not statistically significant, given the small size of our sample, but 
it represents indeed a trend, compatible with the the results present in the literature. More research work is necessary 
to be done to confirm the presence of this interesting effect. 



/'• 


Id 


EO 


EC 


HI 


1.730 (0.010) 


1.800 (0.010) 


H2 


1.740 (0.012) 


1.770 (0.013) 


BP 


1.595 (0.005) 


1.746 (0.018) 


Dl 


1.748 (0.013) 


1.787 (0.014) 


D2 


1.960 (0.023) 


1.974 (0.024) 



TABLE I: This Table shows a sample of the results of the analysis on our group of individuals in the EO (left side of the 
Table) and EC (right side of the Table) conditions. Results show a larger value of \x for EC condition. In parenthesis we report 
the standard deviation. The symbol H indicates a healthy individual, BP the individual with back pain, D the subjects with 
mental depression. 

Before ending this Section, we want to stress that the main result of this paper, namely the surprising complexity 
matching between brain and complex sounds that we shall discuss in Section [Vj does not depend in any way on the 
conjecture we make about the emergence of fi < 2 in the single electrodes. This conjecture would have the effect of 
explaining the findings of Refs. [flEEl]. but, if proved wrong, would not weaken the validity of the main result of this 
paper. 
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EO 


EC 


1.261 


1.888 


1.242 


1.931 


1.294 


1.898 


1.278 


1.877 



TABLE II: This Table shows a sample of the results of the analysis of Ref. [45| on healthy individuals in the EO (left side of 
the Table) and EC (right side of the Table) conditions. The parameter /j, has been obtained according to Eq. (|14|l . Results 
show a larger value of \i for EC condition. 

V. ANALYSIS OF MUSIC COMPOSITION 

The music composition data that we analyze in this Section have been produced b y m eans of a virtual instrument 
(software synthesizer) suitably designed [13, EH to produce physiological effects [4!j[5(3]. The instrument is capable 
of generating abstract sonic textures that are free from overt cultural influences. Moreover, through the use of 
Presets, a record can be kept of all parameters of consequence in the generation of the musical material. This 
feature is important for purposes of correlating the output sonic textures with, for instance, biological data from 
EEG recordings. In this way, researchers can precisely pinpoint areas of interest in biological data for analysis of and 
comparison with the generative parameters of the sound structures. Designed to allow control over the redundancy 
of time-point and pitch/frequency patterns in a hierarchical framework, the virtual instrument sonic textures can be 
gradually morphed between constant states (stable) regimes and chaotic (complex) regimes via controls built into the 
graphical user interface. Moreover, precise measurement and recording of all generative schemes is possible, as well as 
synchronization with the EEG data time series. An important characteristic of the instrument is that vertical sonic 
textures are flexibly configurable with respect to the degree of vertical and horizontal redundancy of pitch- and time- 
space organization within the ongoing sonic flow. Various degrees of complexity can be introduced via Presets prior 
to and during the transmission of audio data to the test subject. For instance, microtonal controls are built into the 
instrument, so that the researcher can regulate the content of the sonic texture at any given point to create a relative 
harmonicity shift of the composite sonic data stream. In other words, there is maximum control over the degree of 
complexity via overlapping, phase-modulated patterning of melodic, harmonic, rhythmic, and texture-density data. 

The instrument produce a vector V(t), with 7 components, namely pitch onset time (in milliseconds), note frequency 
(in Hertz), waveform type, amplitude, articulation, preset number, and oscillator number. All these components are 
filed and assumed to afford information about the signal. We therefore adopt for the music composition the same 
procedure applied before for the EEG signal, building also in this case a time series of MST topological changes. The 
resulting SP for a sample record is plotted in Fig. [3] 

We see from the figure that also in this case the distribution of topological changes produces a non exponential SP, 
namely a stretched exponential with parameters 7 = 0.430 and a = 0.600. As in Section ITVl a further test is needed 
to prove that the process is renewal. The results of the AE are plotted in Fig. [U 

Also in this case the AE supports the renewal assumption for the process. We conclude that music composition 
shares the same properties as the EEG signal analyzed in Section ITVl This suggests that both human brain and music 
composition, sharing the condition /i < 2, are complex systems that, in the absence of the exponential truncation, 
would violate the ergodic condition [5lj . 

VI. A SHORT REVIEW OF THE CM EFFECT 

The authors of Refs. [H, have developed a linear response theory that applies to non-Poisson renewal systems. 
The form of this linear response theory is given by: 

n(i) = <&(*)) = e [ dt' X {t,t')t P {t'), (15) 
Jo 

where £s(t) denotes the signal produced by the system of interest, x(A^') is the linear response function and £p(i) 
is the external perturbation. In the case of the human brain the signal £s(£) is a global property emerging from 
the synchronization of different brain areas, roughly corresponding to the superposition of the signals detected by 
the electrodes fixed on the patient's scalp. In the absence of perturbation this global signal is characterized by the 
non-Poisson renewal events revealed by the method described in Section [TVl To transform these data into we 
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FIG. 3: The SP of MST topological changes obtained from the music composition. The best fitting function is a stretched 
exponential of the form of Eq. @, with 7 = 0.430 and a = 0.600. t a = 100. 
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FIG. 4: The AE on the time series of Fig. [3] The aging effect is compatible with the renewal assumption. 

assign alternate signs to the quiescent time regions between two consecutive events. The complexity of this signal is 
denoted by means of the power index /is < 2, emerging from the analysis of Section ITVl 

If the perturbation function £ p (t) is harmonic the system does not respond [231 ) . This is a clear sign of complexity 
insofar as a non-Poisson renewal signal cannot be interpreted as the superposition of infinitely many independent 
processes. The individual constituents of a set of neurons responsible for any cognitive action are expected to be 
strictly cooperating with all the others' constituents. As a consequence, a harmonic perturbation triggers a cascade 
of different time scales, thereby violating the prescriptions of ordinary stochastic resonance processes [23|, [53[. Thus, 
we are led to make the conjecture that the transmission of information from the perturbing signal to the human 
brain is determined by the interaction between the renewal events of the perturbing signal and the renewal events 
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of the perturbed system. If the perturbing signal does not have any renewal events, as in the case of a harmonic 
perturbation, there is no response to a weak perturbation. 

More recently the authors of Ref. [2(J HH have proved that this conjecture is correct, and that in the case where 
£p(i) is a signal derived from another non-Poisson renewal system with index \xp < 2, the system responds, and 
the intensity of the response is maximum when we use the matching condition /is — /ip[20l|. The authors of Ref. 
[20L l2lj] have denoted this effect with the name of Complexity Matching (CM), and they proved that when /ip < /is 
the perturbed system inherits the perturbation power index. In Section |V] we have proved that music composition is 
actually a complex signal with /ip < 2. Thus, the existence of the CM effect leads us to conjecture that the reason 
why the brain is sensitive to music lies on the fact that both the brain and music are non-Poisson renewal systems 
living in the non-ergodic region. 



VII. CONCLUDING REMARKS 



The literature on complexity is wide, and there are many different proposals to account for complex processes. All 
these proposals share only one essential property, this being the departure from the canonical exponential distribution 
of ordinary statistical physics. The authors of Ref. [4l[ have pointed out that this departure can be realized by 
means of quite different physical processes, either non-homogeneous Poisson processes or homogeneous non-Poisson 
processes. The different physical origin of these two processes is revealed by the AE. 

According to the authors of Refs. [H, [5f| there exists a close connection between self-organized criticality [f|, 
superstatistics (56j and non-extensive thermodynamics [57| . This analysis coincides with the critical view illustrated 
by Jensen [58l |. It is very attractive to conjecture that the brain operates at or near a self-organized critical state 
[55l . [59( 1 . This corresponds to the recent observation [43], [6(| that neuron synchronization is a sort of phase transition 
involving a close cooperation among the elementary constituents of the neuron set. However, in our opinion, the main 
limits of these interesting theories is that they do not pay attention to the important role of renewal events, whose 
objective existence is made compelling by the results of the AE, in both the case of blinking quantum dots [4jJ and 
of neuron synchronization (43j . 

We think that, although these properties have been overlooked by the majority of the researchers working in the 
field of complexity, they deserve more attention, and we hope that this paper may serve the important purpose of 
raising the interest of the investigators in this direction. 

If the importance of these non-Poisson renewal events is not recognized, the results of the analysis of this paper on 
both EEC data and music composition, which are based on a solid method of statistical analysis [41] , are incompre- 
hensible. If, on the contrary, we accept the leading idea that complexity rests on the close cooperation of elementary 
components losing their own identity at the moment of the onset of synchronization, then the fact that both music 
composition and the human brain are non-Poisson renewal processes becomes a natural way of explaining why music 
exerts its influence on the brain. 

This conclusion is so important as to deserve further remarks. First of all, a more appropriate term to denote the 
music composition analyzed in this paper would be complex sound. In fact, the main purpose of the music composition 
utilized in this work is to affect brain complexity rather than generate emotional and/or aesthetic responses, per se. 
Thus, to make the composition of these complex sounds more flexible and more suitable to the purpose of realizing 



an efficient transport of information from acoustic excitation to the human brain, thanks to the CM effect 2^, 



we 



21 



2fl|21| 



have to focus on the realization of a given /ip more than on esthetic purposes. On the basis of the CM effect 
the network of EEC electrodes is expected to inherit the same fip as that of the acoustic excitation if /ip < /is, where 
Us denotes the brain complexity index. We plan to design a real experiment based on recording EEG data from a 
subject listening to music with an exponent /Ltp.We plan also to observe whether the exponent of the EEG changes if 
the listeners hears songs with different exponents. All of this requires much more work, and the solution of technical 
as well as conceptual problems. However, we are convinced that the results of this paper are important enough as to 
trigger further research work for the realization of this important experiment. 



APPENDIX A 



Let us address the problem of establishing the SP ^>sp(t) corresponding to a set of sequences prepared at t = 0. 
It is evident that 

00 „t 

* S p(t) = V/ dt'Mt')f(t-t')exp(-rn). (Al) 
n=o Jo 
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To explain Eq. (|A1[) , let us notice that 'fsp(i) is the probability that the physical generator does not produce any 
further event, after the initial preparation event, until time t. In the natural time scale the probability density of 
not producing an event is exp(-rn). To evaluate Vl'sp(t), we have to multiply exp(-rn) by i/) n (t') and by ^(t — t'). 
This is because exp(-rn) indicates that until time t the physical generator acts n times, thereby making n drawings 
from the distribution density ^(t). Note that ip(t) is also called subordinating function (32J. This function is the 
probability distribution density affording information on the time distance between two consecutive actions of the 
physical generator, not necessarily producing events. With no drawing (n = 0), we get ^sp(t) = &(t), where 'P(t) is 
the SP corresponding to namely, the probability that the physical generator does not act until time t. With n 
drawings we fill a time interval of length T\ + T2 + ■ ■ -r„ = t' < t. The factor ^(t — £') ensures that no drawing occurs 
in between t' and t. The function ip n (t') denotes the probability that n drawings from the distribution ip( T ) occurred, 
the last of which occurred exactly at time £'. Due to the renewal nature of this process we have 

Mt) =^n-l(t)®Mt), (A2) 

where ® indicates time convolution and ipi(t) — ip(t). 

The Laplace transform of 9sp(t), &sp(u), is expressed as a function of the Laplace transform of ip(t), V>(u), as 



follows 



= T : — T - \l - i){u) . (A3) 



1 — ip{u)exp{— r) 



u 



To obtain the previous expression, still no hypothesis on the form of ip(r) has been done. If ^(t) is an exponential, 
the subordination exerts no physical effect on the SP, as its form remains exponential. Let us see the case of non 
exponential subordination function. With straightforward algebra, and assuming r -C 1, Eq. (|A3[) becomes: 



that is the Laplace transform of 



= ATT ^ = 1 — T7T ( A4 ) 

u + r<P(u) 1 — ip(u) 



j t *sp(t) = -r j * S p(t - t')$(t')dt' (A5) 

and <&(£') is a memory kernel. Note that the assumption r -C 1 is not necessary to generate the time convoluted 
structure of Eq. (IA5j) . It is essentially required to make the natural time n compatible with a continuous time 
representation, thereby yielding for the survival probability in the natural time scale the exponential form exp(-rn). 
An even more consequence of r 1 is that, as we shall see hereby, this condition makes very extended the stretched 
exponential regime. 

In the case where the memory kernel $(i) is a delta of Dirac, Eq. (|A5|) makes 'l'sp(t) become an ordinary 
exponential. To generate a stretched exponential we must make a proper choice of the memory kernel $(i), and 
consequently of the subordination function. Let us assign to the memory kernel in the Laplace space the following 
form 

Hu) = x (u + r) 2 -» X ^[r(2- A1 )T^ 1 ]- 1 . (A6) 

and let us assume the parameter \x to be smaller than 2. With this choice the SP becomes, in the Laplace domain, 

*sp(u) = - 1 , ru _ a (A7) 

with 

a = n-l j=( X r)i. (A8) 

If T = 0, we recognize in Eq. (|A7j) the well known Laplace transform of a ML function of order a, see Eq. ([3|). For 
a = 1, the ML function becomes an exponential. 

Note that the parameter T has the important role, as in Ref. [32| . of taking into account that we are working with 
systems of finite, rather than infinite size. Consequently, we must assign to T a finite value. Thus, let us assume that 
r > and that r < 7 < 1. Note that 7 <C 1 is generated by r <C 1. The condition T < 7 is made necessary by the 
request that the truncation of the fat tail of the subordinating function leaves some sign of the system complexity. 
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However, this condition implies a departure from the pure ML relaxation function, and the form of "Jsp on the time 
scale at which we observe the histogram and SP. In the short-time regime, 7 <C u <C 1, it is impossible to neglect the 
first term in the denominator of Eq. (|A7j) . and therefore, considering that r <C u 

^sp(u) = - I 1 _ a , (A9) 

which is the Laplace transform of a stretched exponential. On the contrary, if the condition r « a < 7 applies, the 
first term in the denominator of Eq. (|A7|) can be neglected, and the expression for the SP reads 

which, thanks to the Tauberian theorem, is the Laplace transform of an inverse power law SP. In this condition, the 
fat ML function tail becomes visible. If T < 7, namely, T is moderately smaller than 7, Eq. (|A7|) shows that the 
inverse power law never appears. 

Let us now find analytically the form of the subordination function ip(t). If we adopt the expression for the memory 
kernel of Eq. (|A6|) , we obtain 

»'°4gT°rpr° » ' < A11 » 

Again, two regimes clearly appear depending on the parameters involved in Eq. (|A11[) : if V <C u <C 1, then Eq. (|A11[) 
becomes 

ip(u) ~ 1 - -u"- 1 (A12) 

X 

that is the Laplace transform of an inverse power law distribution density. If we explore the regime u < L, Eq. (|A11[) 
reads 

- 1 T 1 5 ( Al3 ) 

J- T X p2- M 

that is the Laplace transform of an exponential distribution density, implying that the subordination function ij)(t) is 
an inverse power law, truncated exponentially in the long time limit. 

In conclusion, in the case of systems of finite size the fat tail of the subordination function is truncated. In the case 
where r is not very small, and consequently 7 is only moderately small, the stretched exponential regime is not very 
extended and the adoption of a truncated subordination function generates a \E'sp(t), with a distinct inverse power 
law tail. If, r <C 1, and consequently the stretched exponential regime is very extended, a truncated subordination 
function may have the effect of canceling the inverse power law tail of the SP, and the stretched exponential remains 
the only sign of complexity. However, the subordination function has an inverse power law nature, and its power 
index /1 is derived from a through Eq. (|A8[) . thereby yielding Eq. ([4]). 
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